Integer Quantum Hall Effect in Trilayer Graphene 
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The Integer Quantum Hall Effect (IQHE) is a distinctive phase of two-dimensional electronic 
systems subjected to a perpendicular magnetic field. Thus far, the IQHE has been observed in 
semiconductor heterostructures and in mono- and bi-layer graphene. Here we report on the IQHE 
in a new system: trilayer graphene. Experimental data are compared with self-consistent Hartree 
calculations of the Landau levels for the gated trilayer. The plateau structure in the Hall resistivity 
determines the stacking order (ABA versus ABC). We find that the IQHE in ABC trilayer graphene 
is similar to that in the monolayer, except for the absence of a plateau at filling factor v — 2. At 
very low filling factor, the Hall resistance vanishes due to the presence of mixed electron and hole 
carriers induced by disorder. 

PACS numbers: 61.72.Bb, 71.55.Cn 



More than 30 years after its initial discovery in semi- 
conductor two dimensional electron gases (2DEG), the 
Integer Quantum Hall Effect (IQHE) remains one of 
the most fascinating phenomena in condensed matter 
physics. The recent discovery of graphene [l[ boosted 
this research field by providing a new 2D system where 
Dirac-like electronic excitations with Berry's phase ir 
leads to a new form of IQHE 0, H|, with plateaus at 
<?xy = ( n + |) <?e 2 /h, where g is the Landau level degen- 
eracy due to spin and valley degrees of freedom. Soon 
afterward, a third type of IQHE was reported in bilayer 
graphene, where the 27r Berry's phase of charge carriers 
results in a conventional quantization sequence, except 
that the last Hall plateau is missing [|| . As the dynamics 
of charged carriers change every time an extra graphene 
layer in added, it was theoretically anticipated that the 
Landau Level (LL) spectrum of A^-layer graphene sys- 
tems would result in distinctive IQHE features arising 
from an A^7r Berry's phase. In trilayer graphene, the 
zero-energy LL is expected to be 12-fold degenerate so 
that the Hall resistance plateau sequence follows the same 
ladder as in graphene, but the plateau at v = ±2 should 
be missing (see figures QJa and QJb). So far, most of 
the studies dedicated to IQHE in trilayer graphene have 
been restricted to theoretical calculations [5|-[9|. Indeed, 
experimental realizations are difficult since the knowl- 
edge of the exact number of layers and their relative 
stacking order are challenging to ascertain. Making use 
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of both high field magneto-transport and Raman spec- 
troscopy, we clearly identified a trilayer graphene sam- 
ple and report on a fourth type of IQHE in this sys- 
tem. Self-consistent Hartree calculations of the gated 
trilayer Landau levels based on the Slonczewski-Weiss- 
McClure (SWMC) tight binding model [IM1 have been 
performed and favorably compared to the experimental 
data. Thus, a comprehensive knowledge of the under- 
lying electronic properties of trilayer graphene is pre- 
sented, together with an unambiguous determination of 
the stacking order. 

Many few layer graphene flakes were deposited onto a 
d = 280 nm thick thermally grown silicon oxide on sili- 
con substrate using the standard micro-mechanical exfo- 
liation of natural graphite. The flake measured in this 
study was first roughly identified as few-layer graphene 
using optical microscopy. The Raman scattering spec- 
trum was measured at room temperature using a con- 
focal micro Raman scattering set-up using He-Ne laser 
excitation (A = 632.8 nm) with ~ 1 mW optical power 
focused on a 1 [im diameter spot. The 2D band feature 
(also called G' feature) of this sample is shown in figure 
[T}c and appears in the form of a multicomponent feature 
characteristic of multi-layer graphene specimens [l3j , dif- 
ferent from the one observed for a mono-layer graphene 
specimen processed in the same way (also shown in fig- 
ure [T]-c). The experimental IQHE of the sample under 
study is displayed in figure [l}d, together with bi-layer 
and mono-layer graphene fingerprints of other samples. 
These samples have an equivalent carrier density and 
similar mobility (see legend of figured] for details). For 
the sample under consideration, the sequence of the Hall 
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FIG. 1: Identifying the sample under study as a trilayer graphene. a) and b) schematic representations of the IQHE in mono-, 
bi- and tri-layer graphene. The degeneracy of the zeroth Landau Level, equally shared by electrons and holes, is / = g-fold 
degenerate in graphene, / = 2g-fo\d degenerate in bilayer graphene and / = 3g-fold degenerate in trilayer graphene where 
g — 4, leading to a sequence of the Hall resistance plateaus successively shifted by 2h/e 2 Q.. Possible LL degeneracy lifting at 
high field is not represented in these drawings, c) Raman spectrum of trilayer graphene and mono-mayer graphene measured on 
the same substrate, d) Experimental IQHE in tri, bi and mono-layer graphene for equivalent carrier density n = 3.4 x 10 12 cm~ 2 
and quasi-equivalent mobility. The optical image of the trilayer graphene sample is shown in insert. As contact 1 was proven 
defective, a constant current of i — 1 /iA is injected through contacts 2 and 4. The Hall resistance is measured between contacts 
3 and 5. 



resistance plateaus is described by R xy = h/ve 2 where 
v = 6, 10, 14.... Taking into account the results of both 
Raman spectroscopy and IQHE, we thus unambiguously 
identify the sample as trilayer graphene. This assertion 
is further reinforced when comparing the optical contrast 
of the flake with other mono-layer graphene flakes under 
microscope (see the supplementary information). As ex- 
pected, the IQHE in trilayer graphene is indistinguish- 
able from its monolayer counterparts except at very high 
field where the v = 2 quantized Hall resistance plateau 
is absent. 

Before further proceeding with a detailed presenta- 
tion of IQHE in graphene trilayer, let us review some 
important experimental considerations. Once identified, 
the sample was connected to multi-terminal electrodes 
made by successive thermal evaporation of Ti [5 nm] 



and Au [40 nm] through a PMMA mask (see figure 
[TJd, inset). The flake was further patterned in the 
Hall bar geometry using electron beam lithography 
and oxygen plasma etching. The contact resistances 
were estimated to a few hundred ohms each. A gate 
voltage V g between the sample and the substrate was 
used to electrically induce electrons or holes up to 
densities of 6 x 10 12 cm -2 . The sample was annealed 
alternatively in vacuum and in helium gas at a tem- 
perature T=110°C. The induced carrier density as a 
function of the gate voltage follows the simple plane 
capacitor model n = k€q{V s — VcNp)/ed where en is the 
vacuum dielectric permittivity, n = 3.9 is the relative 
permittivity for Si02 and Vcnp = —13.75 V is the gate 
voltage required to reach the Charge Neutrality Point 
(CNP), indicating the presence of a n-type residual 
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FIG. 2: a) Electrical resistance of the sample as a function 
of carrier concentration at T = 4.2 K. T heoretical adjust- 
ment of the experimental data using |l4| together with the 
deduced field-effect mobilities at residual carrier density of 
0.85 x 10 12 cm -2 b) Quantized Hall resistance profile in tri- 
layer graphene for selected charge carrier density far from 
CNP. 



doping estimated to n = 0.85 x 10 12 cm~ 2 . Figure [2]-a 
shows the typical two-probe longitudinal resistance as 
a function of the carrier density at T = 4.2 K, while 
figure HJ-b displays the Hall resistance for selected carrier 
concentrations far from CNP. The lower bound (Hall) 
mobility is estimated to fi c ~ 1300 ± 100 cm 2 .V _1 ■ s _1 
for electrons and fj,h ~ 900 ± 300 cm 2 . V -1 -s _1 for holes. 
The Hall resistance displays well defined plateaus at 
R xy = h/{i>e 2 ) with v = 6, 10, and 14. We notice 
that the Hall resistance slightly overshoots the v = 6 
resistance plateau at very high magnetic field. This 
surprising feature triggered the need for a detailed 
theoretical analysis of the LL spectrum [l5|, [l6| which 
goes beyond the simple model presented earlier in the 
introduction. 

To model the gated trilayer, we employ an SWMC 
tight-bindin g p aramctrization of the local hopping am- 
plitudes [13, Hll and treat the Coulomb interactions via a 
self-consistent Hartree approach. There are two possible 
stacking orders to consider: Bernal (ABA) and rhom- 
bohcdral (ABC), both illustrated in figure |31 For the 



Bernal case, we take 70 = 3000 mcV, 71 = 400 meV, 
72 = — 20mcV, 73 = 300 meV, 74 = 150 meV, and 
75 = 38 meV. In addition, there is an on-site energy 
shift of |A = 18meV for each c-axis neighbor. For the 
rhombohedral case [l2| , the parameter 75 does not enter. 
When a gate voltage is applied to the device, the charge 
carriers will be distributed among the three layers. In- 
tegrating Gauss' law across the layers provides a set of 
three equations which have to be solved self-consistently 
in order to obtain the trilayer potentials (see supplemen- 
tary informations for details). We assume that the sys- 
tem contains stray charges which are located on the top 
layer of the device. These charges are responsible for 
the offset voltage Vqnp necessary to bring the system to 
neutrality. The zero temperature Hall conductivity shifts 
by e 2 /h when the electrochemical potential crosses the 
center of a Landau level. We include Zecman splitting, 
assuming g = 2. Finally, some disorder is introduced in 
the model through a convolution of the density of states 
with a square distribution of half- width Wq = 10 meV, 
estimated from the sample's mobility. 

Figure [3] shows the LL energies and theoretical Hall 
resistance for both ABC and ABA stacking for gate volt- 
age V g = +50 V, along with the experimental results 
for R xy (the full set of IQHE data is reported in the 
supplementary information). For fields up to 40 T, the 
measurements agree fairly well with the theoretical pre- 
dictions for the ABC trilayer, and fail to reproduce the 
theoretical Hall plateau sequence for the ABA trilayer. 
The contrasting plateau sequences for ABC and ABA 
trilayers arise due to the significant differences in their 
respective Landau level structures, as is evident in figure 
[3] Indeed, the rhombohedral stacking order accounts for 
the absence of Hall plateaus at some filling factors, like 
v — 8 and v = 12. A similar effect occurs in monolayer 
graphene, where the plateaus at v = 4, 8, 12, ... are 
missing due to valley degeneracy arising from the inver- 
sion symmetry of the honeycomb lattice. With no bias 
voltage, the ABC trilayer is inversion symmetric, while 
the ABA trilayer is not [17j . The presence of an electric 
field across the graphene layers, due to the gate voltage 
induced charge redistributions [H, [l8| , breaks the lattice 
inversion symmetry. However, this field-induced splitting 
between LLs arising from different valleys is small in the 
ABC case, except for the six levels closest to the Dirac 
energy. Neglecting Zeeman splitting, quantum Hall steps 
of amplitude Aa xy = 2.e 2 /h should be observed for each 
plateau-to-plateau transition. This holds in particular for 
the Bernal type of stacking as the LLs originating from 
valleys K and K' are quite distinct from each other due 
to the absence of inversion symmetry. On the other hand, 
the ABC-stacked LL band-structure is much less affected 
by electrostatic effects. In high enough magnetic field, 
the LLs evolve roughly by bunches of four and, when 
disorder effects are taken into account, lead to quantum 
Hall steps of Aa xy = 4.e 2 /h, as experimentally observed. 
It has been demonstrated that the stacking in graphene 
trilayers can be analyzed making use of micro- Raman 
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FIG. 3: a) Illustration of ABC-stacked trilayer 
graphene (top panel), theoretical Landau level 
structure (middle panel) and quantized Hall resis- 
tance (bottom panel) using the following parame- 
ters : V e = 50 V, Vcnp = -13.75 V, T = 4.2 K, 
Wo = 10 meV and g = 2. In the middle panel, the 
solid and dashed curves indicate the Landau levels 
originating from valleys K and K', respectively. In 
the bottom panel, the experimental Hall resistance 
is displayed as the red curve, superimposed with 
theoretical results as the black line, b) For com- 
parison, the same information for the ABA stacked 
trilayer graphene 
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spectroscopy fl9| . However, the Raman spectral signa- 
ture of stacking is quite subtle and undoubtedly incon- 
clusive for the present experiment (figure [TJc). Although 
it has been found that ABC-stacked graphene trilayer 
flakes are much rarer than their ABA-stackcd counter- 
part using the micro-mechanical cleavage method, the 
theoretical analysis of IQHE allows an unambiguous de- 
termination of the rhombohedral staking order for the 
sample considered here. 

Interestingly, the IQHE fails to be reproduced at 
very low filling factor (high magnetic field and low 
carrier density). To further investigate this issue, we 
analyze the Hall resistance for various charge carrier 
concentrations close to CNP. Focusing at figure QJa, 
we begin the analysis with the Hall resistance for 
Vg = +40 V (n = 4.3 x 10 12 cm~ 2 ), which displays well 
defined quantized Hall plateaus. As the gate voltage 
is decreased, the corresponding set of curves is shifted 
so that the quantized plateaus occur at lower magnetic 
field, as expected from IQHE theory. On the other hand, 



as the Fermi energy is driven closer to the CNP, the low 
field Hall effect is no longer linear reflecting the presence 
of electrons and holes that both contribute to transport. 
In the range —50 V < V g < +20 V, the initial ratio 
between electron/hole density evolves as the magnetic 
field is increased to accommodate the field-induced 
redistribution of quantum states available in the lowest 
LL [2l[. Actually, the electron and hole densities tend 
to equilibrate and consequently the Hall resistance 
vanishes at high field. This effect is a hallmark of the 
disordered 2DEG (22), where the presence of electron 
and hole puddles allows both types of carriers for a 
given Fermi energy close to CNP. Indeed, the inevitable 
presence of charged impurities or lattice imperfections 
introduces long-range disorder and leads to a spatially 
inhomogeneous and fluctuating potential landscape, 
resulting in some local accumulation of charge carriers 
[23T - I271 ] . Alternatively, its consequences on charge trans- 
port can be monitored through temperature-dependent 
measurements. The main frame of figure HJ-b shows the 
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FIG. 4: a) Hall resistance for various back gate voltage in 
the close vicinity of the CNP. Notice that the Hall resistance 
tends to vanish at very high magnetic field as the ratio be- 
tween electrons and holes tends to equilibrate, b) Longitudi- 
nal resistance R2-4 as a function of gate voltage and temper- 
ature. The maximum resistance as a function of temperature 
is shown in the insert. A theoretical fit according to the model 
developed in [2(J is also shown. The theoretical adjustment 
takes into account the increase of carriers density with tem- 
perature and the resulting activated conductance through an 
inhomogeneous system made of electron/hole puddles. 



two-probe longitudinal resistance -R2-4 with varying 
gate voltage and temperature. The resistance increases 
with decreasing temperature in the vicinity of the CNP 
while it remains constant for higher doping. This trend 
is consistent with what has been previously observed 



in disordcrcd-likc samples, for which the highly inho- 
mogeneous carrier density can lead to both metallic 
and activated transport. The temperature evolution of 
the resistance maximum at CNP is given in the inset 
of figure 0Jb. A good agreement with the predictions 
of reference |20| is achieved with fitting parameters 
(i = 1150cm 2 /V-s and A$ = 58meV, where A$ 
stands for the rms amplitude of the fluctuating potential 
landscape. This value is slightly larger than those 
reported in the literature for bilayer graphene [2^, [29[ 
(A$ = 21.5 meV for \x w 2000cm 2 /V-s), as predicted for 
FLG [H. 

To summarize, we report for the first time the observa- 
tion of a new form of IQHE in a gated graphene trilayer. 
The filling factor sequence associated with the quantized 
Hall resistance plateaus is identical to that for graphene, 
but the plateau at v = 2 is missing. The experimental 
data are supported by a theoretical analysis where 
both Bernal and rhombohedral stacking have been 
considered. The main experimental IQHE features are 
reproduced only for the rhombohedral case, emphasizing 
the importance of stacking order in the electronic prop- 
erties of graphene trilayers. Usually, graphene trilayers 
deposited on Si02 show a poor mobility of the order of 
1000cm 2 /V-s ^Mtt^, justifying the need for very high 
magnetic field for IQHE studies. Important progresses 
are expected in samples with higher mobility, e.g. in 
suspended or boron nitride-deposited graphene trilayers. 
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